*, Nobuo
Ueda2, Martin Gühmann3, Luis Alfonso
Yañez-Guerra1, Piotr Słowiński1, Kyle C. A.
Wedgewood1 and Gáspár
Jékely1,4,*: k.jokura @ exeter.ac.uk,
gaspar.jekely @ cos.uni-heidelberg.deNitric oxide (NO) produced by nitric-oxide synthase (NOS) is a key regulator of animal physiology. Here we uncover a function for NO in the integration of UV exposure and the gating of a UV-avoidance circuit. We studied UV/violet avoidance mediated by brain ciliary photoreceptors (cPRCs) in larvae of the annelid Platynereis dumerilii. In the larva, NOS is expressed in interneurons (INNOS) postsynaptic to cPRCs. UV stimulation of cPRCs triggers INNOS activation and NO production. NO signals retrogradely to cPRCs to induce their sustained post-stimulus activation through an unconventional guanylyl cyclase. This late activation inhibits serotonergic ciliomotor neurons to induce downward swimming. In NOS mutants, retrograde signalling, circuit output and UV avoidance are defective. By mathematical modelling, we recapitulate phototransduction and circuit dynamics in wild-type and mutant larvae. Our results reveal how NO-mediated retrograde signalling gates a synaptic circuit and induces short-term memory of UV exposure to orchestrate light-avoidance behaviour.
In nervous systems, synaptic transmission and volume transmission together shape circuit dynamics (Bargmann and Marder, 2013). While synaptic transmission occurs at specialised contact sites, volume transmission is characterised by the delocalised release of diverse diffusive neuromodulators.
Nitric oxide (NO) is one such modulator with unique physical and signalling properties. This free radical synthesized from L-arginine by nitric oxide synthase (NOS) is short-lived and can diffuse across biological membranes (Cudeiro and Rivadulla, 1999; Thomas, 2015). Canonical NO signalling involves the Ca2+/calmodulin-dependent activation of NOS, NO production and diffusion, and the NO-dependent activation of soluble guanylate cyclases (sGC). This can occur cell-autonomously or in other cells leading to cGMP production (Bredt et al., 1990; Hölscher, 1997). Given that NOS activation is calcium dependent and NOS shows neuron-type-specific expression (Aso et al., 2019; Gibbs and Truman, 1998; Mobley et al., 2022; Wildemann and Bicker, 1999), NO action can lead to the activity-dependent modulation of neural circuits at specific sites (Aso et al., 2019; Jacoby et al., 2018; Vielma et al., 2014; Wang et al., 2007) .
In the vertebrate retina, NOS is expressed in amacrine, ganglion and other cells and the actions of NO can be diverse (Cudeiro and Rivadulla, 1999; Jacoby et al., 2018; Wang et al., 2007). For example, defective NO signalling in NOS knockouts leads to a decreased sensitivity of retinal ganglion cells to light stimulation (Wang et al., 2007). In the retina, NO signalling can also involve pathways other than canonical sGC signalling (Jacoby et al., 2018; Tooker et al., 2013; Wei et al., 2012). Due to the complex expression of NOS in vertebrates and the diversity of its functions, it has been challenging to link the neurophysiological effects of NO to signalling mechanisms and behaviour change.
In the Drosophila brain, NO signalling is also involved in tuning circuit activity at diverse and specific sites. In the ellipsoid body of the central complex, NO is involved in visual working memory (Kuntz et al., 2017). NO signalling also tunes the dynamics of associative memories in mushroom body circuits. NOS is expressed in PPL1-γ1pedc mushroom-body neurons where it is involved in shortening memory retention while promoting fast memory updating in response to new experiences (Aso et al., 2019).
In the whole organism context, NO signalling has often been studied in diverse marine invertebrates, where NO can regulate larval settlement and metamorphosis (Leise et al., 2001; Locascio et al., 2022; Song et al., 2021; Ueda et al., 2016; Zhang et al., 2012). While NOS-expressing neurons potentially responsible for these effects have been reported (Bishop and Brandhorst, 2007; Locascio et al., 2022), it has not been possible to link these NO-dependent effects on behaviour or life-cycle transitions to neuronal activity and function.
Overall, we still know little about how NO production relates to stimulus conditions, how it shapes circuit activity at specific neuron types, and how NO-dependent modulation relates to behaviour.
To investigate NO function in neural circuit dynamics and behaviour, here we study larvae of the marine annelid model Platynereis dumerilii (Ozpolat et al., 2021). Platynereis has emerged as a model for systems neuroscience with a toolkit that enables the combination of behavioural analysis and neuronal activity imaging with genetic manipulations. A whole-body synaptic connectome and gene expression atlases are also available and can be integrated with functional approaches in live animals thanks to the cellular-level stereotypy of larvae of the same developmental stage (Ozpolat et al., 2021; Verasztó et al., 2020; Verasztó et al., 2017; Vergara et al., 2021).
Here we uncovered an essential function for NO signalling in larval UV/violet-light avoidance behaviour. In Platynereis larvae, UV/violet avoidance is mediated by brain ciliary photoreceptor cells (cPRCs) and is characterised by downward swimming (Verasztó et al., 2018). The cPRCs express a ciliary-type opsin, c-opsin1 (Arendt et al., 2004) that forms a UV-absorbing bistable photopigment with an absorption maximum around 384 nm (Tsukamoto et al., 2017; Veedin Rajan et al., 2021; Verasztó et al., 2018). Upon UV/violet exposure, the cPRC show a characteristic biphasic calcium response that is c-opsin1 dependent. UV/violet avoidance is also c-opsin1-dependent and is defective in c-opsin1 mutants (Verasztó et al., 2018). Here we show that Platynereis NOS is expressed in interneurons of the cPRC circuit and is required for UV/violet-avoidance. By combining calcium imaging across the fully-mapped cPRC circuit (Verasztó et al., 2018) with genetic perturbations and mathematical modelling, we describe how NO tunes circuit dynamics through non-synaptic retrograde signalling to cPRCs. This delayed neuroendocrine feedback integrates UV/violet exposure to induce a short-term memory manifested in altered circuit activity and an aversive behavioural response.
Nitric oxide synthase is expressed in interneurons of the UV-avoidance circuit
We identified a single nitric oxide synthase (NOS) gene in the Platynereis dumerilii genome and transcriptome data. Phylogenetic analysis of NOS proteins indicate that Platynereis NOS belongs to an orthology group of bilaterian NOS sequences (Figure 1—figure supplement 1). To characterise the expression pattern of NOS we used in situ hybridization chain reaction (HCR) and transient transgenesis. In two- and three-day-old larvae, we detected NOS expression in four cells (two of them weakly expressing) in the apical organ region (Figure 1D and Figure 1—figure supplement 2). NOS was also expressed in the region of the visual eyes (adult eyes) and the pigmented eyespots (Figure 1—figure supplement 2). The four apical organ cells, but not the eyes, were also labelled with a NOS reporter construct driving palmitoylated-tdTomato (Figure 1E). This reporter also revealed the axonal projections of these central NOS-expressing neurons. The position and morphology of the four NOS+ cells allowed us to identify the same four cells as four interneurons (INNOS) in our three-day-old whole-body Platynereis volume EM data (Verasztó et al., 2020; Williams et al., 2017) (Figure 1B,C). In the synaptic connectome, the INNOS cells are postsynaptic to the UV-sensory cPRCs and presynaptic to the INRGW interneurons, which are also cPRC targets (Figure 1C,F). The projections of INNOS cells are segregated into input (dendritic) and output (axonal) compartments, with cPRC inputs occurring on the dendritic part in the dense neurosecretory plexus. INNOS output synapses form in the more ventral projection region (Figure 1—figure supplement 3). INRGW neurons synapse on the head serotonergic ciliomotor neurons (Ser-h1), which synapse on the prototroch ciliary band and the cholinergic MC ciliomotor neurons (Figure 1C,F) (Verasztó et al., 2017).
Figure 1. Identification of NOS-expressing interneurons (INNOS) within the cPRC circuit. (A) Scanning electron microscopy image of a three-day-old Platynereis larva. (B, C) Volume rendering of the neuron types (cPRC, INNOS, INRGW, Ser-h1 and MC) in the cPRC circuit reconstructed from a whole-body transmission electron microscopy volume of a three-day-old larva. Neurite skeletons are shown with cell-body positions represented by spheres. Projections of all neurons in the body are shown in grey to highlight the neuropils. The outline of the yolk is also indicated in grey. In B, nuclei positions of the prototroch head ciliary band are shown as grey spheres. (D) Expression of the NOS gene detected by in situ HCR (magenta) in a two-day-old larva (anterior view). Antibody staining for acetylated α-tubulin (acTub: green) highlights cPRC cilia and the neuropil. (E) Expression of a NOS reporter (NOSp::palmi-3xHA: magenta) labelled with an anti-HA antibody in a two-day-old larva (anterior view). Antibody staining for acetylated α-tubulin (acTub: green) highlights cPRC cilia and the neuropil. (F) Synaptic wiring diagram of the cPRC circuit. Hexagons represent cell groups, with the number of cells per group shown in square brackets. Arrows represent the summed number of synaptic contacts between cell groups. Arrow thickness is proportional to the number of synapses.
Nitric oxide is produced during UV/violet stimulation of the cPRCs
The expression of NOS in the INNOS interneurons in the cPRC circuit suggests that NO signalling may be involved in UV/violet-avoidance. To test this, first we asked whether NO is produced during UV/violet stimulation of the larvae. We injected the fluorescent NO-reporter DAF-FM into zygotes and imaged two-day-old larvae while exposing the region of cPRC cilia to 405 nm violet light. To mark cell outlines, we coinjected mRNA encoding a red-fluorescent reporter (RGECO), allowing the identification of the cPRCs. Following light stimulation in the region of the ramified cilia of the cPRCs, we detected an increase in DAF-FM fluorescence in the anterior neurosecretory neuropil, the region of INNOS projections. This increase did not occur in larvae where we illuminated a control area (Figure 2).
Figure 2. NO produced by UV/violet stimulation to cPRCs. (A) DAF-FM fluorescence in the region of the neurosecretory neuropil. The white line indicates the outline of the larva. The dashed line corresponds to the area where fluorescence was quantified. Circles indicate the location of cPRC and control stimulation. cPRCs are marked by thin lines. (B) Changes in DAF-FM fluorescence before and during 405 nm light stimulation. (C) Changes in DAF-FM fluorescence over time during 405 nm stimulation of the cPRCs (green) or a control area (ctr: gray). The purple box indicates the duration of 405 nm stimulation. Individual traces normalized (ΔF/F0) are shown as thin lines. Thick lines show the mean value with 0.95 confidence intervals. N = 9 larvae for control and 11 for cPRC stimulation.
Nitric oxide signalling mediates UV-avoidance behaviour
We next tested whether NO signalling regulates UV/violet avoidance. To achieve this, we generated two Platynereis NOS knockout lines with the CRISPR/Cas9 method. We recovered two deletions (NOSΔ11/Δ11 and Δ23/Δ23), both frame-shift mutations leading to an early stop codon and thus likely representing null alleles (Figure 3—figure supplement 1A). We could establish a homozygous line for both mutations indicating that NOS is not an essential gene in Platynereis. To quantify UV avoidance, we recorded the trajectories of freely swimming wild type and mutant larvae in vertical columns, illuminated laterally from two opposite sides with 395 nm UV light (Figure 3A and Figure 3—figure supplement 1B). As previously shown, wild-type larvae swim downward following non-directional UV/violet light stimulation (Verasztó et al., 2018). In contrast, both two- and three-day-old homozygous NOS-mutant larvae showed a strongly diminished UV-avoidance response (Figure 3A, B and Figure 3—figure supplement 1B,C). This phenotype is similar to the defective UV-avoidance of c-opsin1 mutant larvae (Verasztó et al., 2018) and reveals a requirement for NOS in UV-avoidance behaviour. Wild type but not mutant larvae also showed an increase in swimming speed under UV light that may be due to downward swimming trajectories (swimming in the direction of gravity) (Figure 3B and Figure 3—figure supplement 1C). We also tested directional phototaxis, by exposing larvae to 480 nm directional collimated light from the top of the column. Three-day-old but not two-day-old NOS-mutant larvae also showed reduced phototactic behaviour, suggesting a function for NOS in the visual eyes that mediate three-day-old phototaxis (Randel et al., 2014) (Figure 3D and Figure 3—figure supplement 1G).
To distinguish between an acute and developmental function of NOS in light responses, we next tested larvae exposed to the NOS inhibitor L-NAME. Larvae incubated for 5 min in 0.1 mM or 1 mM L-NAME showed a dose-dependent inhibition of UV avoidance. In contrast, phototaxis was not affected (Figure 3C, E). Overall, our results indicate an acute requirement for NOS signalling in UV-avoidance and a possible indirect, developmental role in the visual system, reminiscent of the function of NO signalling in Drosophila eye development (Gibbs and Truman, 1998).
Figure 3. NOS is required for UV avoidance in Platynereis larvae. (A) Swimming trajectories of wild type (WT, n=32) and NOS mutant (NOSΔ11/Δ11, n=26 and NOSΔ23/Δ23, n=47) three-day-old larvae. All trajectories start at 0 x and y position and time 0 corresponding to 10 sec after the onset of 395 nm stimulation from the side. (B) Vertical position of batches of wild type and mutant three-day-old larvae over time under 395 nm UV stimulation. The starting position of each larval trajectory was set to 0. (C) Vertical position of batches of control and L-NAME-treated (0.1 and 1 mM) three-day-old larvae over time under 395 nm UV stimulation. The starting position of each larval trajectory was set to 0. (D) Vertical displacement in 30 sec bins of wild type and mutant (NOSΔ11 and NOSΔ23) three-day-old larvae stimulated with 395 nm light from the side, 488 nm light from the top and 395 nm light from the top. (E) Vertical displacement in 30 sec bins of control and L-NAME-treated (0.1 and 1 mM) three-day-old larvae stimulated with 395 nm light from the side, 488 nm light from the top and 395 nm light from the top.
NO retrograde signalling tunes cPRC responses to UV/violet stimulation
To investigate how NO signalling alters the dynamics of the cPRC circuit, we carried out calcium imaging experiments. We ubiquitously expressed the calcium sensor GCaMP6s in larvae and imaged calcium signals during 405 nm stimulation of the cPRCs. As we have shown previously, a 20-sec local stimulation of cPRC cilia led to a transient increase in cPRC calcium levels, followed by a transient decrease (Verasztó et al., 2018). After ~20-sec, calcium levels in cPRCs were raising again, reaching higher levels than at the start of the stimulus – a response that may involve depolarisation (Figure 4A, B). This activation phase occurs after the 20 sec stimulation period and is likely due to a delayed neuroendocrine feedback (Verasztó et al., 2018). To determine whether NO mediates such a feedback, we repeated the experiment in NOS-mutant larvae. While we detected the initial activation phase followed by inhibition, in homozygous NOS-mutants for both CRISPR alleles this was not followed by delayed activation. Instead, calcium levels dropped to a low steady-state level (Figure 4A, B). We thus identified a requirement for NO signalling in the late-phase activation of cPRCs.
Two unconventional guanylyl cyclases are expressed in the cPRCs
We aimed next to identify the NO receptor in the cPRCs. NO generally acts via soluble guanylate cyclases (sGC), belonging to the guanylate cyclase family with a CYC domain (PFAM domain: PF00211). NO binding to the heme group of sGC leads to increased cyclic guanosine monophosphate (cGMP) production. Analysis of sGCs in Platynereis indicated that these genes are not expressed in any of the cells of the cPRC circuit (Verasztó et al., 2017). Recently, Moroz and coworkers reported an atypical but widely conserved family of guanylyl cyclases with a NIT (nitrite/nitrate sensing) domain (PF08376) (NIT-GC) as potential mediators of NO signalling (Moroz et al., 2020). To identify NIT-GCs in Platynereis, we searched transcriptome resources and retrieved 15 potential NIT-GC homologs (Figure 4—figure supplement 1 and 2). To analyse the relationship of these sequences to metazoan NIT-GCs, we retrieved protein sequences with a CYC domain from the transcriptome and genome databases of 45 metazoan and 2 choanoflagellate species. We carried out cluster analysis and did phylogenetic reconstruction on a group of membrane-bound guanylyl cyclases with sGCs as an outgroup. In agreement with Moroz et al. (Moroz et al., 2020), we found a group of GCs with NIT domains with representatives in placozoans, cnidarians, some ecdysozoans, echinoderms, and lophotrochozoans. The 15 Platynereis sequences belonged to several deeply diverged clades in the phylogenetic tree (Figure 4—figure supplement 1 and 2).
To characterise the expression of NIT-GCs, we used previously published spatially mapped single-cell transcriptome data (Achim et al., 2015; Williams et al., 2017). Among the 15 NIT-GCs, two showed high and specific expression in the cPRCs and one was expressed in the INNOS cells (Figure 4—figure supplement 2). In the single-cell data, we could identify the cPRCs by the specific expression of c-opsin1 and the pedal-peptide2 neuropeptide precursor (MLD proneuropeptide), previously described cPRC markers (Arendt et al., 2004; Williams et al., 2017) (Figure 4—figure supplement 3A). The INNOS cells were identified by NOS expression and spatial mapping in the brain (Achim et al., 2015). We decided to focus on two NIT-GCs expressed in the cPRCs and with a full-length sequence, NIT-GC1 and NIT-GC2. To confirm the single-cell data, we first carried out in situ hybridisation chain reaction (HCR) with probes for NIT-GC1 and NIT-GC2 mRNA. Both genes were specifically expression in the four cPRCs, as confirmed by co-labeling with an acetylated α-tubulin antibody and with an HCR probe against pedal peptide 2/MLD proneuropeptide (Figure 4C, D and Figure 4—figure supplement 3A-C). To analyse the subcellular localisation of NIT-GC1 and NIT-GC2 at the protein level, we raised and affinity-purified polyclonal antibodies against a specific peptide sequence from both proteins. In immunostainings, we found that NIT-GC1 was localise to the region corresponding to the axonal projections of the cPRCs in the anterior neurosecretory plexus (Figure 4E). Co-immunostaining with the rabbit NIT-GC1 antibody and a custom rat antibody raised against Platynereis NOS revealed the localisation of both proteins in dots in the neurosecretory neuropil (Figure 4—figure supplement 3D). In contrast, NIT-GC2 specifically labelled the ramified sensory cilia of the cPRCs (Figure 4F). These different subcellular localisations suggest that the two NIT-GCs are involved in different intracellular signalling processes in the ciliary and axonal regions of the cPRCs.
NIT-GC1 produces cGMP in an NO-dependent manner
To further characterise these atypical guanylyl cyclases, we focused on NIT-GC1 and carried out in vitro experiments. In bacteria, NIT domains are thought to regulate cellular functions in response to intra- or extracellular nitrate and nitrite. NIT-GC1 has a NIT domain and a highly conserved cyclase domain that is expected to catalyse cGMP synthesis (Figure 4G). The NIT domain may render NIT-GC1 dependent on NO signals. To test this, we co-expressed the cGMP indicator Green cGull (Matsuda et al., 2016) and NIT-GC1 in cultured COS-7 (monkey kidney) cells, a cell line with minimal endogenous sGC activity (Matsuda et al., 2016). For balanced expression, we used a single plasmid with the two open-reading frames separated by the 2A self-cleaving peptide (Figure 4H). Application of the NO donor SNAP lead to increased Green cGull fluorescence, an effect we did not observe when cells were exposed to DMSO or when Green cGull was expressed alone (Figure 4I-K). To test whether this effect is dependent on the NIT domain, we also tested a deletion construct of NIT-GC1 lacking the NIT domain (Figure 4G). Cells expressing this construct and Green cGull did not show an increased fluorescence of the cGMP reporter when exposed to SNAP (Figure 4L). These results indicate that NIT-GC1 is able to catalyse cGMP production in an NO-dependent manner and this function requires the NIT domain. These results establish NIT-GC1 as a biochemical sensor of NO or its derivatives.
NIT-GC1 is required for NO-mediated retrograde signalling to cPRCs during the UV response
To test the in vivo function of NIT-GC1 and NIT-GC2 in cPRC responses, we combined calcium imaging with morpholino-mediated knockdowns. We used two translation-blocking morpholinos for each NIT-GC gene and tested knockdown efficiency by immunostaining injected animals with the NIT-GC1 and NIT-GC2 antibodies (Figure 4—figure supplement 3E,F). For both genes, the morpholinos led to a strong reduction in the respective antibody signal, confirming efficient knockdown and antibody specificity.
In NIT-GC1 morphant larvae, the delayed activation of cPRCs following 405 nm stimulation did not occur (Figure 4M). This phenotype is similar to the phenotype of NOS mutants suggesting that NIT-GC1 acts as the NO sensor in cPRCs to drive their delayed activation. This could occur via increased cGMP production and the opening of a cyclic-nucleotide-gated ion channel (CNG) specific to cPRCs (Tosches et al., 2014). NIT-GC2 morphant larvae, in contrast, showed a step-up increase in calcium following light stimulation (Figure 4N). The calcium signal decayed during stimulation and was off after light off. These data support an essential role for ciliary-localised NIT-GC2 in suppressing cPRC calcium following its transient rise at stimulus onset. Overall, these knockdown experiments revealed different signalling mechanisms for the two NIT-GCs that may be due to their different subcellular localisations.
Figure 4. NOS and two NIT-GCs shape calcium signals during cPRC UV/violet response. (A, B) GCaMP6s signals in cPRCs in wild type and NOS mutant (A, NOSΔ11/Δ11, B, NOSΔ23/Δ23) larvae during 405 nm light stimulation. (C, D) In situ HCR for (C) NIT-GC1 and (D) NIT-GC2 (magenta) in three-day-old Platynereis larvae. Larvae were co-stained with an antibody against acetylated α-tubulin to label cPRC cilia and the neuropil (green). (E, F) Immunostaining for (E) NIT-GC1 and (F) NIT-GC2 (magenta), co-stained for acetylated α-tubulin (green). (G) The domain structure of Platynereis NIT-GC1 and the truncated NIT-GC1ΔNIT protein lacking the NIT domain. A predicted transmembrane region (TM) is shown in grey. (H) Schematic of the cell-based assay to detect cGMP production following the addition of an NO donor SNAP or DMSO as control. (I-L) Green cGull fluorescence over time for the four conditions tested. Individual responses and their mean with 0.95 confidence interval are shown (n > 6 cells). Intensities are normalized (ΔF/F0). The indicated chemicals were added at 2 min after the start of imaging (grey bars). (M, N) GCaMP6s signals in cPRCs in (M) NIT-GC1 and (N) NIT-GC2 morphant larvae during 405 nm light stimulation. Individual responses and their mean with 0.95 confidence interval are shown.
NO signalling shapes the dynamics of the cPRC circuit
To investigate how NO and NIT-GC signalling influence the dynamics of the cPRC circuit, we imaged calcium signals from postsynaptic neurons in wild type, mutant and morphant larvae. We were able to image the activity of all neurons in the cPRC circuit (INNOS, INRGW, Ser-h1 and MC). The MC cell was identified based on its position and intrinsic activity (Verasztó et al., 2017). To unambiguously identify all other cells from which we recorded calcium signals, we developed an on-slide immunostaining method (Figure 5—figure supplement 1A). We used the cell-specific antibody markers against RYamide (INNOS) (Figure 5—figure supplement 1B-D), RGWamide (INRGW) and serotonin (Ser-h1) (Conzelmann et al., 2011) to immunostain agar-embedded larvae following calcium imaging. Based on the position of the nuclei, we could correlate live and fixed samples at a single-cell precision (Figure 5A,B). Due to the stereotypy of the larvae, we could also identify neurons based on their position and calcium activity in activity-correlation maps (Figure 5C).
We first quantified the responses of the INNOS and INRGW interneurons during 405 nm stimulation of the cPRCs. In both wild type and NOS-mutant larvae, INNOS cells showed an increase in calcium during stimulation (Figure 5D). In contrast, the INNOS response was flat or slightly negative in NIT-GC2 morphant larvae (Figure 5E) revealing an essential role for NIT-GC2-mediated cPRC suppression in INNOS activation. INRGW cells were initially inhibited during cPRC stimulation, followed by a delayed activation paralleling the second activation phase of cPRCs. This late INRGW response was lacking in NOS-mutants (Figure 5F). In NIT-GC2 morphants, INRGW cells showed a transient increase in calcium that decayed after light off and a delayed activation was not present (Figure 5G).
Next, we imaged calcium signals from Ser-h1 and MC neurons in wild type and NOS mutant larvae. Ser-h1 cells showed an activation profile that correlated with cPRC activity, including a reduction in calcium during stimulation followed by rebound, a response that was defective in NOS mutants (Figure 5I). MC cells showed sustained activation, including a late-phase that was lacking in NOS mutants (Figure 5I). These data suggest that during 405 nm stimulation the Ser-h1 cells are inhibited and MC cells are activated, and this regulation is NO-dependent. This pattern is expected to inhibit ciliary activity in the prototroch but not in the other ciliary bands, triggering NO-dependent downward swimming.
Figure 5. NOS- and NIT-GC2-dependent dynamics of the cPRC circuit. (A, B) GCaMP6s imaging from cPRCs and INNOS cells (left panels) followed by on-slide immunostaining for (A) RYamide to label INNOS and (B) RGWamide+serotonin to label INRGW and Ser-h1 (red). Nuclei are stained with DAPI (cyan). Asterisks indicate cPRC nuclei. Numbers mark the same cells in the GCaMP and immunostaining images matched by position. (C) Correlation map of neuronal activity of the cPRCs, INNOS, INRGW, Ser-h1 and MC neurons. (D) GCaMP6s fluorescence in INNOS cells in wild type (WT) and NOSΔ11/Δ11 mutant larvae during 405 nm stimulation of the cPRC cilia. (E) GCaMP6s fluorescence in INNOS cells in NIT-GC2 morphant larvae during 405 nm stimulation. (F) GCaMP6s fluorescence in INRGW cells in wild type and NOSΔ23/Δ23 mutant larvae during 405 nm stimulation. (G) GCaMP6s fluorescence in INRGW cells in NIT-GC2 morphant larvae during 405 nm stimulation. (H, I) GCaMP6s fluorescence in (H) Ser-h1 cells and (I) the MC cell in wild type and NOSΔ11/Δ11 mutant larvae during 405 nm stimulation.
Mathematical modelling of cPRC-circuit dynamics
To further analyse the dynamics of responses to UV light and formally describe cPRC phototransduction, synaptic connections and NO retrograde signalling, we developed a mixed cellular-circuit-level mathematical model. We used our Ca2+ imagining data of cPRC, INNOS and INRGW cells collected in wild type, NOS knockout, and NIT-GC2 morphant larvae to formulate assumptions that are the basis of the proposed model.
We model a direct c-opsin1-dependent response to UV leading to the initial rise in Ca2+ levels. The mechanism of this signal is not known but may include Gα and Gβ signalling, CNGAα channel activation or other pathways. We further assume that in cPRC cells, Ca2+ levels increase with cGMP. We model a decrease in cGMP and Ca2+ levels due to a NIT-GC2-dependent pathway (Figure 6A). The effect of NO in cPRC cells is captured by another term that describes a NIT-GC1 and NO-dependent increase in cGMP. Based on the synaptic connectome, we infer a feedforward coupling between cPRC cells and INNOS and INRGW cells. The proposed synaptic coupling is UV-dependent, decays linearly and is suppressed in NIT-GC2 morphants (Figure 6A).
In the INNOS cells, we assume that excitatory synaptic input from cPRC (or rebound from tonic inhibition) leads to a rise in Ca2+ leading to NO production via NOS. This input is NIT-GC2-dependent. In the INRGW cells, we model a decrease in Ca2+ based on an inhibitory UV-dependent synaptic signal from the cPRCs. We further assume the existence of direct excitatory coupling between cPRC Ca2+ levels and INRGW Ca2+ to account for the late effects of cPRC on INRGW. The UV-dependent and Ca2+-dependent signals may be mediated by different transmitters released by cPRCs during distinct phases of their activation cycle. In addition, we assume a direct inhibitory coupling between INNOS Ca2+ and INRGW Ca2+. Finally, in all variables, we assumed a linear decay and constant production to set a steady state (Figure 6A). Since the aim of the model is to capture the normalised fluorescence data, the model is nondimensionalised. UV stimulation is modelled as a square pulse. To find model parameters producing an output fitting the experimental data we employed a global optimisation method known as a genetic algorithm (GA).
Our model only includes a minimal set of parameters and assumptions of interactions that are required to describe the dynamics of cPRC, INNOS and INRGW in wild type and loss-of-function conditions (Figure 6A, B). The model highlights that cPRC phototransduction employs distinct pathways that operate on different time scales and differentially influence cPRC Ca2+ levels. The coupling between cPRCs and interneurons also requires different signalling mechanisms, either through distinct neurotransmitters or receptors expressed in the different cells.
To identify possible molecular pathways, we analysed single-cell transcriptome data for cPRC, INNOS and INRGW identified by spatial mapping (Achim et al., 2015) and unique marker genes (Williams et al., 2017). In cPRCs, we found transporters and synthetic enzymes indicating cholinergic, glutaminergic, glycinergic, GABAergic and adrenergic neurotransmission (Figure 6C). For each neurotransmitter, we also found receptor-encoding genes expressed in INNOS and INRGW (Figure 6C). The two types of interneurons often expressed different subunits or types of these receptors, indicating differential signalling (Figure 6C).
Based on these data, our experimental results and the mathematical model, we assembled a minimal phototransduction and circuit diagram. The components for which experimental data are available are shown in bold. Other potential molecular players are also indicated (Figure 6D).
Figure 6. Mathematical modelling and signalling mechanisms of the cPRC circuit. (A) Diagram of the mathematical model with componenets, interactions, parameters and equations used to model Ca2+ dynamics. (B) Average calcium traces and modelled traces fitted to the data. (C) Dot plot of genes (columns) expressed in three types of cells (rows) in the cPRC circuit using single cell RNA-Seq. The size of the dots is expressed in proportion to the percentage of cells expressing that gene relative to all cells. The colours represent the normal logarithm of the number of transcripts in the cells expressing the gene. (D) Schematic diagram of the signalling pathway of the cPRC circuit, focusing on the NO feedback.
Our work revealed an essential role for NO-mediated signalling in driving UV/violet avoidance behaviour in larval Platynereis. NO, produced by postsynaptic INNOS interneurons, signals retrogradely to presynaptic cPRCs via NIT-GC1 leading to delayed and sustained cPRC activation. This late-phase activation of the photoreceptors drives circuit output through projection interneurons and ciliomotor neurons. In the cPRC circuit, synaptic connectivity alone is thus not sufficient to account for circuit dynamics and behavioural change, as documented in other circuits (Bargmann and Marder, 2013; Imambocus et al., 2022).
Localised NO signalling in the neurosecretory plexus
NO is a free radical with a millisecond-to-second half life and thus a limited signalling range. In neuronal signalling, NOS is often localized to neurites (Kuntz et al., 2017) at close proximity to sGC at synapses (Burette et al., 2002; Garthwaite, 2015). In the cPRC circuit, NOS is localised to the dendritic compartment of INNOS cells where also cPRC postsynaptic sites occur. NIT-GC1, the target of NO signalling is localised to cPRC projections. NO-mediated retrograde signalling thus likely occurs in the neurosecretory plexus where NOS- and NIT-GC1-containing projections are in close proximity and where we detected NO production following UV stimulation. In contrast, INNOS to INRGW synapses occur outside the neurosecretory plexus in the ventral brain neuropil.
Compartmentalised signalling also occurs in peptidergic modulatory systems and can enable selective network activity during specific behaviours. In the UV-avoidance circuit of Drosophila larvae, a peptidergic hub neuron Dp7 links UV-sensory neurons (v’td2) and motor circuits. Dp7 expresses an an insulin-like peptide Ipl7 that is required for acute and sustained UV avoidance. Ilp7 signalling occurs at a functionally and morphologically distinct dendritic compartment of Dp7, segregated from other sensory-motor pathways involving Dp7 (Imambocus et al., 2022).
Functional diversity of NIT-GCs
NO signalling is commonly mediated by sGCs. We identified 12 sGCs in Platynereis, but none of these is expressed in the cPRC based on the available scRNAseq data. Instead, we identified an unconventional cPRC-expressed NIT-domain containing GC, NIT-GC1 as the mediator of NO retrograde signalling. In an in vitro assay, we could show that NIT-GC1 can produce cGMP following the addition of an NO donor and that this activity requires the NIT domain.
NIT domains were first identified in bacteria and animal NIT-GCs have only recently been reported (Moroz et al., 2020; Shu et al., 2003). Bacterial NIT domains regulate cellular functions in response to changes in extracellular and intracellular nitrate and/or nitrite concentrations (Camargo et al., 2007). NO is readily converted to nitrate and nitrite (Garthwaite, 2015; Möller et al., 2019; Santos et al., 2011) and these molecules accumulate in placozoans and cnidarians in cells and tissues with high NOS activity (Moroz et al., 2020, 2004). NIT domains in NIT-GCs may also sense nitrate and nitrite, as in bacteria, a possibility we cannot rule out based on our cellular assays with NIT-GC1. If different NIT-GCs have different sensitivities to NO, nitrite and nitrate, then a range of activation timings may be possible due to the different half-lives of these molecules (Lundberg et al., 2011).
NIT-GC1 and NIT-GC2 showed specific cellular co-expression but very different subcellular localisation and function. In Platynereis, we identified 15 NIT-GCs, suggesting a wide range of functions. Differences in subcellular localisation and biochemical function thus seem to also contribute to the diversity of NIT-GC functions in addition to differences in expression.
Mechanism of phototransduction and neurotransmission in the cPRC circuit
Based on our data herein and previous work we can now propose a more detailed model of cPRC phototransduction and neurotransmission. The cPRCs have high basal Ca2+ and respond to UV/violet light dependent on c-opsin1 (Verasztó et al., 2018). c-opsin1 forms a bistable photopigment and signals through Gi/oα and Gβγ (Tsukamoto et al., 2017; Tsukamoto and Kubo, 2023; Veedin Rajan et al., 2021). In heterologous systems, the Gβγ subunits released following c-opsin1 activation open GIRK channels inducing K+ efflux (hyperpolarisation) (Tsukamoto et al., 2017; Tsukamoto and Kubo, 2023) and close voltage-gated Ca2+ channels, thereby reducing intracellular calcium levels (Tsukamoto and Kubo, 2023). A GIRK channel is also expressed in the cPRCs (Figure 6A). The pathway for the first rapid cPRC activation phase following UV/violet stimulus is not known but may involve the activation of a CNGα channel expressed in the cPRCs (Tosches et al., 2014). For the second inhibitory phase of phototransduction, we identified a key requirement for ciliary-localised NIT-GC2. The mechanisms may involve c-opsin1-dependent inhibition of tonic NIT-GC2 activity and the reduction of ciliary cGMP.
NIT-GC2-dependent signalling is required for the feedforward activation of the INNOS cells through an unknown transmitter. The activation of INNOS leads to NOS activation and NO release, potentially through a canonical Ca2+-calmodulin pathway. Our model predicts feedforward inhibiton from INNOS to INRGW, possibly mediated by glycin transmission. NO released by INNOS neurites activates NIT-GC1 in cPRC projections that could lead to cGMP production and the opening CNGAα. This late-phase cPRC activation is not directly dependent on the UV/violet signal and can happen post-stimulus. The late-phase cPRC activation leads to INRGW activation via an unknown transmitter that is likely different from the one acting on INNOS.
In addition, the cPRC circuit expresses several neuropeptides and their receptors, suggesting further neuromodulatory signals. For example, INRGW express the proneuropeptide RGWamide and cPRCs and INNOS express its receptor, suggesting retrograde peptidergic signalling in the circuit.
We have recently shown that the cPRC circuit also mediate responses to hydrostatic pressure via the same motor system involving Ser-h1 and the prototroch ciliary band. Increased pressure induces cPRC activation and a circuit output that is inverted relative to UV-induced activation. Consequently, ciliary beating increases and the larvae swim upwards. The effect of pressure on cilia requires synaptic transmission by the serotonergic Ser-h1 neurons (Bezares-Calderón et al., 2023).
Pressure-induced Ca2+ transients in cPRCs lack an inhibitory phase and late activation and resemble UV responses in NIT-GC2 morphants (Bezares-Calderón et al., 2023). This indicates that the differentiation of sensory cues by the multisensory cPRCs occurs already at the level of sensory signal transduction. The different signalling pathways then likely result in the differential release of transmitters and modulators that are decoded by the postsynaptic interneuron circuit. The complex transmitter phenotype of cPRCs could underlie such differential signal processing.
Nitric oxide confers short-term memory to circuit activity
Retrograde signalling by NO from INNOS to cPRC leads to the sustained activation of cPRCs and postsynaptic neurons even after the end of stimulation. This activated state is maintained for several tens of second. NO signalling thus induces a transient circuit state or short-term memory trace in the Platynereis larval brain. Because of the short life time of NO, this molecule may be well suited to encode transient memory traces (Kuntz et al., 2017).
In the ellipsoid body of the Drosophila central brain, NO signalling has a similar function. Here, NOS is specifically expressed in the R3 ring neurons and is required for the short-term (~4 sec) visual memory of objects (Kuntz et al., 2017). NO is produced in the axons of R3 neurons and acts directly on sGC in the same axons. This autocrine signal leads to a CNG-dependent temporary increase in calcium levels, carrying the working memory trace (Kuntz et al., 2017).
In the Platynereis circuit, our mathematical model indicates that the onset, magnitude and duration of the NO-dependent signal depends on the intensity and duration of the UV/violet stimulus. This suggests that the NO-dependent memory trace also encodes the magnitude and duration of the stimulus.
During UV avoidance in the planarian Schmidtea mediterranea, neuropeptide signalling has a similar integratory function. Planarians exposed to UV light for >30 sec remain active for extended periods (several minutes) post-stimulation (Bray et al., 2023). If neuropeptide signalling is defective, animals can still respond to UV light, but do not maintain a latent memory state and do not display increased post-stimulus activity.
The possible mechanism of UV-induced downward swimming
UV/violet light induces downward swimming head down (Verasztó et al., 2018). Since stimulus direction is not relevant, swimming direction must be determined by gravity.
Connectome reconstruction and whole-body cell annotation in the three-day-old Platynereis larva has not revealed any balancer organ to sense orientation in the gravity field (Bezares-Calderón et al., 2019; Verasztó et al., 2020). It follows that body orientation must be determined by physical parameters including the buoyancy, centre of mass and shape of the larva as well as differential ciliary activity. Except for ciliary activity, all these parameters are likely invariant during the UV response. Our data suggest the UV-dependent inhibition of prototroch ciliary beating via Ser-h1 inhibition and MC activation (Verasztó et al., 2017).
We hypothesise that the differential beating of prototroch versus trunk cilia causes a head-up or head-down orientation due to physics alone. During the UV response, prototroch cilia beat slower than trunk cilia, resulting in a head-down stable state (‘rear-wheel drive’). In contrast, during the pressure response prototroch cilia beat faster than trunk cilia (Bezares-Calderón et al., 2023), leading to a head-up orientation (‘front-wheel drive’). Testing this hypothesis will require biophysical experiments and mathematical modelling.
UV-avoidance circuits of extraocular photoreceptors
Animals evolved distinct photosensory systems coupled to non-overlapping circuits and guiding unique behavioural responses. These sensory systems can employ different opsin molecules and be tuned to different wavelength of light. The avoidance of noxious UV/violet light is often mediated by extraocular photoreceptors and their circuits, distinct from the pigmented visual eyes. These two types of systems co-exist in Platynereis where the pigmented eyes and eyespots guide phototaxis behaviour with a maximum sensitivity to cyan light (~500 nm) (Randel et al., 2014; Verasztó et al., 2018). Planarians also have pigmented visual eyes mediating phototaxis to cyan light and peripheral extraocular photoreceptors mediating UV avoidance behaviour (Shettigar et al., 2021). Drosophila larvae have cerebral eyes called Bolwig’s organs involved in phototaxis (Kane et al., 2013) and several types of extraocular UV-sensory cells that tile the body wall (Imambocus et al., 2022).
One common feature of UV-avoidance circuits is their sustained post-stimulus activation following UV exposure. This can involve peptidergic signals as in planaria and maggots (Bray et al., 2023; Imambocus et al., 2022) or NO as in Platynereis. Volume transmission is well suited to integrate light exposure and maintain an internal state following noxious stimulation. The amount of modulator released can scale with stimulus intensity or duration and maintain an altered circuit state due to the slower decay of the diffusive signals relative to synaptic transmission.
| Reagent type (species) | Designation | Source or reference | Identifiers | Additional information |
|---|---|---|---|---|
| Strain (Platynereis dumerilii) | NOSΔ11/Δ11 knockout | This paper | Knockout generated by CRISPR/Cas-9-induced gene editing | |
| Strain (Platynereis dumerilii) | NOSΔ23/Δ23 knockout | This paper | Knockout generated by CRISPR/Cas-9-induced gene editing | |
| Strain (Platynereis dumerilii) | NOSΔ11/Δ23 knockout | This paper | Knockout generated by CRISPR/Cas-9-induced gene editing | |
| Cell line (Cercopithecus aethiops) | COS-7 cell | RRID:CVCL_0224 | Angio-proteomie (CAT no. cAP-0203)?? | |
| Biological sample (Platynereis dumerilii) | Wild type Tübingen strain | Other | NCBITaxon:6359 | Jékely lab strain (Tübingen, Exeter) |
| Gene (Platynereis dumerilii) | NOS | This paper | GenBank_Acc#: | |
| Gene (Platynereis dumerilii) | NIT-GC1 | This paper | GenBank_Acc#: | |
| Gene (Platynereis dumerilii) | NIT-GC2 | This paper | GenBank_Acc#: | |
| NOS: Nitric_Oxide_Synthase | To amplify Promoter & Regulatory region | Fwd | NOSProm2ndF0.6BamHI | AGGGATCCCCCAATGCTTTAGCAGTCAGAGGAG |
| NOS: Nitric_Oxide_Synthase | To amplify Promoter & Regulatory region | Rev | GeR1ASCI | AAGGCGCGCCCCACCACCACCTTTGATATCCATGATGCTCACTTCGC |
| NOS: Nitric_Oxide_Synthase | Mutation Check on Exon3 | Fwd | Exon3 Sequence F-27bp | GGTTCATTGGTTTCGATAACATTGCGG |
| NOS: Nitric_Oxide_Synthase | Mutation Check on Exon3 | Rev | Exon3 Sequence R-27bp | CAGAGTCGATCAGTCTGCATATCTCCA |
| NOS: Nitric_Oxide_Synthase | Sequencing primer for Exon3 mutation check PCR product | Fwd | Exon3 Sequence F-2 | GGTGCTCTCCCGGGTACACAA |
| RNA | sgRNA | 5’-TAGGGCAATACTGGCTCCACTC-3’ | ||
| RNA | sgRNA | 5’-AAACGAGTGGAGCCAGTATTGC-3’ | ||
| RNA | pUC57-T7-RPP2-hSpCas9- HA-2XNLS-GFP | plasmid | Bezares-Calderón et al., 2018 | |
| Antibody | Monoclonal Anti-Tubulin, Acetylated antibody | Sigma-Aldrich | Cat#:T6793, RRID:AB_477585 | |
| Antibody | HA-Tag (C29F4), Rabbit mAb | Cell Signaling Technology | Cat#:3724P | |
| NIT-GC1 polyclonal antibody | CYWLLGRKERRPKRRL-amide | This paper | rabbit | Altabioscience |
| NIT-GC2 polyclonal antibody | CTEGSTKEGKKEGQ-amide | This paper | rabbit | Altabioscience |
| NOS polyclonal antibody | CKPSYELQDPWKTYIWRKD-amide | This paper | Rat | Altabioscience |
| Antibody | RYamide neuropeptide antibody | CRY-amide | rabbit | Conzelmann and Jékely, 2012 |
| Antibody | RGWamide neuropeptide antibody | CGW-amide | rabbit | Conzelmann and Jékely, 2012 |
| Antibody | F(ab’)2-Goat anti-Rabbit IgG (H+L) Cross-Adsorbed Secondary Antibody, Alexa Fluor™ 546 | Invitrogen | Catalog # A-11071 | |
| Antibody | Goat anti-Rat IgG (H+L) Highly Cross-Adsorbed Secondary Antibody, Alexa Fluor™ Plus 594 | Invitrogen | Catalog # A48264 | |
| Antibody | F(ab’)2-Goat anti-Mouse IgG (H+L) Cross-Adsorbed Secondary Antibody, Alexa Fluor™ 647 | Invitrogen | Catalog # A-21237 | |
| Recombinant DNA reagent | NIT-GC1 full | comp411593_c0_seq1_309_F | GGTTGAATAATGACAAGCAAGGAGA | |
| Recombinant DNA reagent | NIT-GC1 full | comp411593_c0_seq1_2717_R | GTGCTATCATTTCCAGGTAAATACCC | |
| Recombinant DNA reagent | NIT-GC2 full | Contig2280_66_F | AATATCTAGCGAAGGAGAACACCTCTCTTC | |
| Recombinant DNA reagent | NIT-GC2 full | Contig2280_2763_R | ATGGCCAGTAATAAACCATCAGTGGTTC | |
| Recombinant DNA reagent | pcDNA3.1(+) vector | Invitrogen | Catalog Number: V79020 | |
| Recombinant DNA reagent | Inverse PCR for the insert region of pcDNA3.1(+) | pcDNA3.1(+)_inv_NheI_fwd | CGTTTAAACTTAAGCTTGGTACCGAG | |
| Recombinant DNA reagent | Inverse PCR for the insert region of pcDNA3.1(+) | pcDNA3.1(+)_inv_NheI_rev | CCAGCTTGGGTCTCCCTATAGT | |
| Recombinant DNA reagent | kozak-NITGC2-T2A-Green cGull | NITGC-T2A-fwd | TATAGGGAGACCCAAGCTGGGCCACCATGACCCAGATG | |
| Recombinant DNA reagent | kozak-NITGC2-T2A-Green cGull | NITGC-T2A-rev1 | GCATGTTAGAAGACTTCCTCTGCCCTCATAATCAAACCCCTCTCT | |
| Recombinant DNA reagent | kozak-NITGC2-T2A-Green cGull | NITGC-T2A-rev2 | AGGGCCGGGATTCTCCTCCACGTCACCGCATGTTAGAAGACTTCC | |
| Recombinant DNA reagent | kozak-NITGC2-T2A-Green cGull | NITGC-T2A-rev3 | TGCTCACCATAGGGCCGGGATTCTCCTC | |
| Recombinant DNA reagent | kozak-NITGC2-T2A-Green cGull | cGull-fwd | TCCCGGCCCTATGGTGAGCAAGGGCGAG | |
| Recombinant DNA reagent | kozak-NITGC2-T2A-Green cGull | cGull-rev | ACCAAGCTTAAGTTTAAACGTTACTTGTACAGCTCGTCCATG | |
| Recombinant DNA reagent | Green cGull-T2A-NITGC2 & NITGC2-T2A-Green cGull | NITGC2_seq_743_fwd | AGCCATCTACGAGTGGTA | |
| Recombinant DNA reagent | Green cGull-T2A-NITGC2 & NITGC2-T2A-Green cGull | cGull_seq_385_rev | TGCCCTTCAGCTCGATG | |
| Recombinant DNA reagent | Green cGull-T2A-NITGC2 & NITGC2-T2A-Green cGull | NITGC2_seq_743_rev | TGACTGACGAACCCTCC | |
| Recombinant DNA reagent | Green cGull-T2A-NITGC2 & NITGC2-T2A-Green cGull | NITGC2_seq_394_fwd | AGATATCTTGAAACGGACGA | |
| Recombinant DNA reagent | kozak-NIT1(seq)-T2A-Green cGull | 2A-cGull_invF_L1 | TGACGTGGAGGAGAATCCCGGCCCTATGGTGAGCAAGGGCGAGGAGCTGT | |
| Recombinant DNA reagent | kozak-NIT1(seq)-T2A-Green cGull | 2A-cGull_invF_L2 | GCAGAGGAAGTCTTCTAACATGCGGTGACGTGGAGGAGAATCCCGGCCCT | |
| Recombinant DNA reagent | kozak-NIT1(seq)-T2A-Green cGull | 2A-cGull_invR_L | TCCTTGCTTGTCATGGTGGCCCAGCTTGGGTCTCCCTATAGTGAGTCGTA | |
| Recombinant DNA reagent | kozak-NIT1(seq)-T2A-Green cGull | NIT1-2A_F_L | GAGACCCAAGCTGGGCCACCATGACAAGCAAGGAGATGCCTGTACTCATG | |
| Recombinant DNA reagent | kozak-NIT1(seq)-T2A-Green cGull | NIT1-2A_R_L | TGTTAGAAGACTTCCTCTGCCCTCTATGACTTTTTCTATGCTTTCTTCGG | |
| Recombinant DNA reagent | kozak-NIT1(seq)-T2A-Green cGull | NIT1seq_remo_invF2 | AGCGTGGAGGTGGGCCTAGACGAAAGAGCTGAAAA | |
| Recombinant DNA reagent | kozak-NIT1(seq)-T2A-Green cGull | NIT1seq_remo_invR2 | AGCTCTTTCGTCTAGGCCCACCTCCACGCTGAATA | |
| Recombinant DNA reagent | kozak-NIT1(seq)-T2A-Green cGull | NIT1seq_remo_F2 | GCCGGTCTTGTCGATCAGGATGATCTGGAC | |
| Recombinant DNA reagent | kozak-NIT1(seq)-T2A-Green cGull | NIT1seq_remo_R2 | GTCCAGATCATCCTGATCGACAAGACCGGC | |
| Recombinant DNA reagent | pUC57-NOSp::Palmi-3xHA-tdTomato (plasmid) | This paper | Promoter construct: injected at 250 ng/μl | |
| Recombinant DNA reagent | pUC57-T7-RPP2-tdTomato-P2A-GCaMP6 (plasmid) | This paper | Used for generating tdTomato-P2A-GCaMP6s mRNA | |
| Plasmid | Green cGull | Addgene | Plasmid #86867 | |
| HCR | NOS | Integrated DNA Technologies | ||
| HCR | NIT-GC1 | Integrated DNA Technologies | ||
| HCR | NIT-GC2 | Integrated DNA Technologies | ||
| HCR | RYa-pNP (GenBank accession: JF811330.1) | Integrated DNA Technologies | ||
| HCR | c-opsin1 (GenBank accession: AY692353.1) | Integrated DNA Technologies | ||
| HCR | MLD/pedal2-pNP (GenBank accession: KF515945.1) | Integrated DNA Technologies | ||
| HCR | CNGAα (GenBank accession: KM199644.1) | Integrated DNA Technologies | ||
| fluorescently labeled hairpins | B2-647 | Molecular Technologies | ||
| fluorescently labeled hairpins | B3-546 | Molecular Technologies | ||
| morpholino | NIT-GC1 MO1 | Gene-Tools, LLC | TGCTTGTCATTATTCAACCAGCAAA | |
| morpholino | NIT-GC1 MO2 | Gene-Tools, LLC | TTCAATTAAACCCTCCAGGTTGCTG | |
| morpholino | NIT-GC2 MO1 | Gene-Tools, LLC | AAATGAAGAGAGGTGTTCTCCTTCG | |
| morpholino | NIT-GC2 MO2 | Gene-Tools, LLC | ATATTCATTATGTGAAGAACTTCCA | |
| plasmid for mRNA synthase (GCaMP6s) | BamHI-T7::RPP2(5UTR)-GCaMP6s(AscI-AgeI)-polyA_KpnI_c1 | Plasmid | Bezares-Calderón et al., 2018 | |
| plasmid for mRNA synthase (RGECO1a) | PUC57-T7-PduRPP2(5UTR)-jRGECO1a | Plasmid | Bezares-Calderón et al., 2018 | |
| Chemical compound, drug | SNAP | Sigma-Aldrich | Cat#:M9020 | 500 μM |
| Chemical compound, drug | L-NAME | Sigma-Aldrich | Cat#:M9021 | 501 μM |
| Commercial assay or kit | Phusion Human Specimen Direct PCR Kit | Thermofisher | ||
| Commercial assay or kit | mMESSAGE mMACHINE Sp6 kit | Thermofisher | ||
| Software, algorithm | Golden Gate TAL | Addgene 1000000024 | ||
| Software, algorithm | Effector Kit 2.0, Fiji perl and Fiji scripts for tracking | PMID: 22743772, https://github.com/JekelyLab/Veraszto_et_al_2018 | RRID:SCR_002285 | 0000d2a |
| Commercial assay or kit | QuickExtract | Epicentre,US | Cat#:QE09050 | |
| Commercial assay or kit | MEGAshortscript T7 Transcription Kit | Ambion, ThermoFisher Scientific | Cat#:AM1354 | |
| Commercial assay or kit | mMESSAGE mMACHINE T7 ULTRA Transcription Kit | Ambion, ThermoFisher Scientific | Cat#:AM1345 | |
| Commercial assay or kit | MEGAclear Transcription Clean-Up Kit | Ambion, ThermoFisher Scientific | Cat#:AM1908 | |
| Software, algorithm | Fiji | NIH | RRID:SCR_002285 | |
| Software, algorithm | R Project for Statistical Computing | R Foundation | RRID:SCR_001905 | |
| Software, algorithm | Imaris Version 8.0.0 | Bitplane, UK. | RRID:SCR_007370 | |
| Software, algorithm | CATMAID | DOI:10.1093/bioinformatics/btp266 | RRID:SCR_006278 | |
| Software, algorithm | PhyML | DOI:10.1093/sysbio/syq010 | RRID:SCR_014629 | |
| Software, algorithm | Gblocks | DOI:10.1080/10635150701472100 | RRID:SCR_015945 |
CRISPR-Cas9 Design and Microinjection
Before designing the small guide RNA (sgRNA) for the sgRNA:Cas9 nuclease, splice sites and polymorphic sites in our laboratory culture were identified to avoid them. The sgRNA targeted the third exon of Platynereis dumerilii NOS (target site: 5’-GGGCAATACTGGCTCCACTC-3’). The sgRNA was assembled from two annealed oligonucleotides (5’-TAGGGCAATACTGGCTCCACTC-3’, 5’-AAACGAGTGGAGCCAGTATTGC-3’) forming overhangs for cloning into a BsaI site of the plasmid pDR27456 (Hwang et al. 2013)(42250, Addgene), which contains next to the BsaI site a tracrRNA sequence. The plasmid was then used to PCR amplify DNA (primers: T7, 5’-AAAAGCACCGACTCGGTGCC-3’) for synthesizing the sgRNA. The DNA was purified with the QIAquick PCR Purification Kit (Qiagen). From the DNA, the sgRNA was synthesized with the MEGAshortscript Kit (Thermo Fisher Scientific) and was purified with the MEGAclear Kit (Thermo Fisher Scientific). Cas9-mRNA was transcribed, capped, and polyA-tailed with the mMessage mMachine Kit and the Poly(A) Tailing Kit (both Thermo Fisher Scientific) from a plasmid (pUC57-T7-RPP2-Cas9) containing the Cas9 ORF fused to 169 base pair 5’ UTR from the Platynereis dumerilii 60S acidic ribosomal protein P2. The sgRNA (18 ng/ml) and the Cas9-mRNA (180 ng/µl) were coinjected into fertilized eggs of Platynereis dumerilii wild-type parents according to an established injection procedure (Conzelmann et al., 2013). The eggs were kept at 18°C for 45 min before injection and were injected at 14.5°C. The injected individuals were kept at 18°C for 5 to 8 days in 6-well- plates (Nunc multidish no. 150239, Thermo Scientific) and then cultured at 22°C until sexual maturity. The mature worms were crossed to wild-type worms and the progeny was genotyped, resulting in two founder lines, which were bred to homozygosity.
NOS sequencing and genotyping
For genotyping of the NOS locus, genomic DNA was isolated from single larvae, groups of 6-20 larvae, or from the tails of adult worms. The DNA was amplified by PCR (primers: 5’-GGTTCATTGGTTTCGATAACATTGCGG-3’, 5’-CAGAGTCGATCAGTCTGCATATCTCCA-3’) with the dilution protocol of the Phusion Human Specimen Direct PCR Kit (Thermo Scientific). The PCR product was sequenced directly with a nested sequencing primer (5’-GGTGCTCTCCCGGGTACACAA-3’). A mixture of wild-type and deletion alleles in a sample gave double peaks in the sequencing chromatograms, with the relative height of the double peaks reflecting the relative allele ratio in the sample.
Vertical column setup for measuring photoresponses
Photoresponses of larvae of different ages were assayed in a vertical Plexiglas column (31 mm x 10 mm x 160 mm water height). The column was illuminated from top with light from a monochromator (Polychrome II, Till Photonics). The monochromator was controlled by AxioVision 4.8.2.0 (Carl Zeiss MicroImaging GmbH) via analog voltage. The light passed a collimator lens (LAG-65.0-53.0-C with MgF2 Coating, CVI Melles Griot) before entering the column. The column was illuminated from both sides with light-emitting diodes (LEDs). The LEDs on each side were grouped into two strips. One strip contained UV (395 nm) LEDs (SMB1W-395, Roithner Lasertechnik) and the other infrared (810 nm) LEDs (SMB1W-810NR-I, Roithner Lasertechnik). The UV LEDs were run at 4 V to stimulate the larvae in the column from the side. The infrared LEDs were run at 8 V (overvoltage) to illuminate the larvae for the camera (DMK 22BUC03, The Imaging Source), which recorded videos at 15 frames per second and was controlled by IC Capture (The Imaging Source).
Comparing behavior of wildtype and NOS-knockout 3-day-old larvae
To compare the behavior of wildtype and NOS-knockout larvae at 3 days in the vertical column, the larvae were mixed and left in the dark for 5 min. The larvae were treated with NOS inhibitors for pharmacology. The NOS inhibitors were L-NAME. The larvae were treated with different concentrations in adjacent columns. The concentrations for the NOS inhibitors were control, 1 mM, 0.1 mM. The larvae were recorded for 1 min in the dark followed by exposure to collimated cyan (480 nm) light from the top of the column for 2 min, then 2 min darkness, and finally collimated UV (395 nm) light from the top of the column for 2 min. Stimulus light was provided by the monochromator (Polychrome II, Till Photonics). Scripts are available at https://github.com/JekelyLab/NOS.
NOS Identification and Phylogenetic Analysis
To identify NOS, we obtained a “seed” database of oxygenase domain in Pfam database, PF02898. From these sequences, we produced a Hidden Markov Model (HMM) and used this to mine the 47 metazoan species, 2 choanoflagellate species and 2 filasterea species investigated. HMM models were run in HMMR3 with an e-value of 1e−15. We ran CD-Hit (Fu et al., 2012) to eliminate redundant sequences (at a 80% threshold). We aligned the sequences with MAFFT version 7, with the iterative refinement method E-INS-i. Alignments were trimmed with TrimAl in gappy-out mode (Capella-Gutierrez et al., 2009). To calculate maximum-likelihood trees, we used IQ-tree2 with the LG+G4 model. To calculate branch support, we ran 1,000 replicates with the aLRT-SH-like and aBayes methods (Minh et al., 2020). The sequences used for the phylogenetic analysis are available in Supplementary file 1, the trimmed alignment is available in Supplementary File 2 and the pre-trimmed data in Supplementary file 3.
NIT-GC Identification and Phylogenetic Analysis
To identify NIT-GCs, we obtained a “seed” database of Adenylate and Guanylate cyclase catalytic domain in Pfam database, PF00211. From these sequences, we produced a Hidden Markov Model (HMM) and used this to mine the 45 metazoan species, 2 choanoflagellate species and 2 filasterea species investigated. HMM models were run in HMMR3 with an e-value of 1e−15. We ran CD-Hit (Fu et al., 2012) to eliminate redundant sequences (at a 80% threshold). To identify clusters, we used the convex-clustering option with 100 jack-knife replicates. The NIT-GCs are extremely well conserved in membrane-bound guanylate cyclases and form an easily recognizable cluster. To analyze the phylogeny of NIT-GCs, the cluster containing these GCs together with membrane-bound guanylate cyclases were parsed and used for tree building. We aligned the sequences with MAFFT version 7, with the iterative refinement method E-INS-i. Alignments were trimmed with TrimAl in gappy-out mode (Capella-Gutierrez et al., 2009). To calculate maximum-likelihood trees, we used IQ-tree2 with the LG+G4 model. To calculate branch support, we ran 1,000 replicates with the aLRT-SH-like and aBayes methods (Minh et al., 2020). The sequences used for the phylogenetic analysis are available in Supplementary file 4, the trimmed alignment is available in Supplementary file 5 and the pre-trimmed data in Supplementary file 6.
Single-cell analysis
We used Achim et al. for the single-cell data (Achim et al., 2015). In Williams et al., they used 107 cells as neurons by removing duplicates from Achim et al. single-cell data, so we used those cells (Williams et al., 2017). Since the raw data were read count data, we normalized them to TPM using Python. After that we converted them to log10. From the sum of the expression levels in 107 cells for each gene, We calculated the percentage in each cell. For each cell, we identified them with marker genes. After created the data in Python, plotted it using R dot plots. RPKM calculates the total number of reads per million bp and then divides by the length of each gene, so it is not possible to compare between samples. Instead, TPM first divides by the length of the gene and then divides by the total number of reads per million bases, which allows for more accurate comparisons between samples. In this case we wanted to compare between samples, so we used TPM. The total TPM of each gene between the samples was used to calculate the percentage of expressed genes. The total TPM values for each gene and the percentage of expressed genes are available in Supplementary file 7.
In situ HCR
Larvae were fixed and treated with Proteinase K, according to the conventional WMISH protocol (Tessmar-Raible et al., 2005), with fixation in 4% paraformaldehyde/ PTW (PBS with 0.05% Tween20) for 2 hr at room temperature, and Proteinase K treatment in 100 µg/ml Proteinase K/ PTW for 3 min (Tessmar-Raible et al., 2005). Specifically, for the HCR protocol, samples were processed in 1.5 ml tubes. Probe hybridization buffer, probe wash buffer, amplification buffer, and fluorescent HCR hairpins were purchased from Molecular Instruments (Los Angeles, USA). Hairpins associated with the b2 initiator sequence were labeled with Alexa Fluor 647, and the hairpins associated with the b3 initiator sequence were labeled with Alexa Fluor 546. To design probes for HCR, we used custom software (Kuehn et al., 2021) to create 20 DNA oligo probe pairs specific to P. dumerilii NOS, NIT-GC1, NIT-GC2, RYa-pNP (GenBank accession: JF811330.1), and MLD/pedal 2-pNP (GenBank accession: KF515945.1). The NOS, NIT-GC1 and NIT-GC2 probes were designed to be associated with the b2 initiator sequence, while the RYa-pNP and MLD/pedal 2-pNP probes were designed to be associated with the b3 initiator sequence. For the detection stage, samples were pre-hybridized in 200 µl of probe hybridization buffer for 1 hr at 37°C, and then incubated in 250 µl hybridization buffer containing probe oligos (4 pmol/ml) overnight at 37°C. To remove excess probe, samples were washed 4× with 1 ml hybridization wash buffer for 15 min at 37°C, and subsequently 2× in 1 ml 5× SSCT (5× SSC with 0.1% Tween20) for 5 min at room temperature. For the amplification stage, samples were pre-incubated with 100 µl of amplification buffer for 30 min, room temperature, and then incubated with 150 µl amplification buffer containing fluorescently labeled hairpins (40nM concentration (2ul of 3uM stock in 150ul amplification buffer, snap-cooled as described; (Choi et al., 2018)) overnight in the dark at 25°C. To remove excess hairpins, samples were washed in 1 ml 5× SSCT at room temperature, twice for 5 min, twice for 30 min, and once for 5 min. During the first 30 min wash, samples were counterstained with DAPI (Cat. #40043, Biotium, USA).
Immunohistochemistry
Whole-mount immunostaining of 2 day old Platynereis larvae fixed with 4% paraformaldehyde were carried out using primary antibodies raised against NIT-GC1, NIT-GC2, NOS, RYamide neuropeptide, RGWamide neuropeptide in rabbit, plus a commercial antibody raised against acetylated tubulin in mouse (Sigma T7451). The synthetic peptides contained an N-terminal Cys that was used for coupling during purification. Antibodies were affinity purified from sera as previously described (Conzelmann and Jékely, 2012). Immunostainings were carried out as previously described (Conzelmann and Jékely, 2012). The NOS promoter (fragment sizes: 12 Kb) was amplified and cloned upstream of 3xHA- Palmi-tdTomato. Larvae injected with promoter constructs (ca. 250 ng/ml) were analysed for reporter expression at 3 days post fertilization using an AxioImager Z.1 fluorescence wide-field microscope (Carl Zeiss GmbH, Jena) and immediately fixed for immunostainings. The protocol followed for immunostaining of HA-tagged reporters was recently described (Verasztó et al., 2017). Specimens were imaged with a LSM 780 NLO or LSM 880 with Airysan Confocal Microscope (Zeiss, Jena).
Calcium imaging
For calcium imaging, 49–55 hpf larvae were used. Experiments were performed at room temperature and larvae were immobilised by being embedded in 2.5% agarose filtered artificial seawater between a slide and coverslip spaced with adhesive tape. GCaMP6s mRNA (1 mg/ml) was injected into zygotes as described previously (Randel et al., 2014). Larvae were imaged on a Zeiss LSM 880 with Airyscan (with a C-Apochromat 63X/1.2 Corr - water) with a frame rate of 1.88 frame/sec and an image size of 512 x 512 pixels. The larvae were stimulated in a region of interest (a circle with 50 pixel diameter) with 405 nm lasers controlled by the Bleaching mode. The imaging laser had a similar intensity than the stimulus laser but covered an area that was 10 times larger than the stimulus ROI.
Cell culture experiment
Green cGull was used for the cGMP assay (Matsuda et al., 2016). A full-length Pdum-NIT-GC1 and -NIT-GC2 coding sequences were amplified by PCR starting from a Platynereis dumerilii cDNA library and cloned into the pcDNA3.1(+) vector using the T2A self-cleaving sequence. Cos-7 cells with low expression of endogenous soluble guanylate cyclase were used as cultured cells for gene expression. This cell line was purchased from Angio-proteomie (CAT no. cAP-0203). The Cos-7 cells were maintained at 37 °C in 35mm dishes (Nunc™ Glass Bottom Dishes) containing 3 mL of DMEM, high glucose glutamax medium (Thermo; Cat. No. 10566016) supplemented with 10% fetal bovine serum (Thermo; Cat. No. 10082147). Upon reaching confluency of approximately 85%, we transfected the cells with the plasmid containing Green cGull-T2A-NITGC1. Transfections were carried out with 150 ng of each plasmid and 0.3 μl of the transfection Lipofectamine 3000 Reagent (invitrogen; Cat. No. L3000001). Two days post-transfection, we removed the culture medium and substituted it for fresh DMEM-medium. For single-wavelength imaging experiments, cells in 35-mm dishes were washed twice and imaged in modified Ringer’s buffer (140 mM NaCl, 3.5 mM KCl, 0.5mM NaH2PO4 , 0.5mM S-3 MgSO4 , 1.5 mM CaCl2 , 10 mM HEPES, 2 mM NaHCO3 and 5 mM glucose). Dishes were mounted on a stage heated at 37 °C and imaging was performed using an inverted microscope (LSM880, Zeiss) equipped with an oil-immersion objective lens (UApo/340, 40, NA = 0.17). Images were acquired using a xenon lamp, 460–495 nm excitation filter, 505-nm dichroic mirror and 510– 550-nm emission filter (Zeiss). S-Nitroso-N-acetyl-D, L-penicillamine (SNAP) was purchased from Sigma-Aldrich (St. Louis, MO, USA) . The exposure time of the EM-CCD camera was controlled by the ZEN software (Zeiss). Images were acquired every 15 s for 10 min and stimulation was initiated 2 min after starting image acquisition. Imaging data analysis was performed using ImageJ (National Institutes of Health, Bethesda, MD, USA).
The mathematical model for the wild type data comprises six dynamical equations:
\[\begin{align*} \frac{d C_P(t)}{dt} &= 1 - \delta^{C_P} C_P(t) + K_{UV}^{C_P} UV(t) C_O(t) + K_{G}^{C_P} G(t) , \\ \frac{dC_O(t)}{dt} &= 1 - \delta^{C_O} C_O(t) - K_{UV}^{C_O} UV(t) C_O(t), \\ \frac{dG(t)}{dt} &= 1 - \delta^{G} G(t) + K_{GC2}^G \left(1 - UV(t)\right) G(t) + K_{GC1}^G N(t) , \\ \frac{dS(t)}{dt} &= UV(t) - \frac{S(t)}{\tau_S}, \\ \frac{dC_N(t)}{dt} &= 1 - \delta^{C_N} C_N(t) + K_{S}^{C_N} S(t), \\ \frac{dN(t)}{dt} &= 1 - \delta^{N} N(t) + K_{C_N}^{N} C_N(t), \\ \frac{dC_R(t)}{dt} &= 1 - \delta^{C_R} C_R(t) - K_{S}^{C_R} S(t) + K_{C_P}^{C_R} C_P(t) - K_{C_N}^{C_R} C_N(t) C_R(t), \\ UV(t) &= \begin{cases} a, & t_{UV\: start} \leq t \leq t_{UV\: end}\\ 0, & \text{otherwise} \end{cases} . \end{align*}\]Data were analysed in ImageJ and R. All figures were assembled in R with the cowplot and patchwork packages. All scripts are available at https://github.com/JekelyLab/Joukra_et_al_NOS (commit ). >>>>>>> 9932f067da10604ec00c679a5b10f91decddfa29
This work was funded by the Wellcome Trust (214337/Z/18/Z). This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No 101020792). KJ has been supported by a JSPS Overseas Research Fellowship, LAYG by a BBSRC Discovery fellowship (BB/W010305/1), PS by a Wellcome Trust Institutional Strategic Support Award (204909/Z/16/Z), and KCAW by the EPSRC Hub for Quantitive Modelling in Healthcare (EP/T017856/1).
Figure 1 – figure supplement 1 Phylogenetic tree of NOS by maximum likelihood. Tree robustness was tested with 1000 replicates of ultrafast bootstrap with the aLRT-SH-like and aBayes methods.
Figure 1 – figure supplement 2 Expression analysis of the NOS gene (magenta) using in situ HCR using the larva at three-day-old in dorsal (A) and anterior view (B). (B) Insert: showing two NOS cells (INNOS_dl and INNOS_vl) on left side.
Figure 1 – figure supplement 3 ssTEM reconstruction of the INNOS (A), INRGW (B) and cPRC (C), anterior view. NS plexus; neurosecretory plexus
Figure 3 – figure supplement 1 (A) Top: The domain and Exon/Intron structure of NOS. Bottom: Close-up region showing the genomic locus of NOS gene and the wild-type sequence (WT) targeted by CRISPR/Cas9. The generated mutants (NOSΔ11, NOSΔ23) are also shown. Pink indicates target sites. Gray shows PAM sequences, red shows stop codons. (B) Overlaid trajectories for WT (n=37) and NOS mutant (NOSΔ11/Δ11, n=18 and NOSΔ23/Δ23, n=8) at two-day-old larvae. 0 sec as the starting point. After 10 sec, UV (395 nm) stimulation from the side. (C) The temporal changes in the vertical position of the WT and mutant two-days-old larvae before and after UV stimulation are shown. The starting points of each larval trajectory are set to 0. After UV stimulation is indicated by purple squares. (D, E) The temporal changes in the distance traveled of the WT and mutant in two-day (D) and three-day-old (E) larvae before and after UV stimulation are shown. (F) The temporal changes in the distance traveled of larvae treated with NOS inhibitors, L-NAME in three-day larvae before and after UV stimulation are shown. (G) Vertical swimming in wild-type (WT) and mutant (NOSΔ11 and NOSΔ23) larvae at two-day-old stimulated with UV (395 nm) light from side, blue (488 nm) light from top and UV (395 nm) light from top. The data are shown in 30 s bins.
Figure 4 – figure supplement 1 Cluster analysis of guanylate and adenylate cyclases. Connections are based on blast similarities < 1e-16 as shown on the upper right. Animal groups are colour coded as shown on the upper left. NIT-GCs, NIT domain containing guanylate cyclases; membrane-bound GCs, Membrane-bound guanylyl cyclases; sGCs, soluble guanylate cyclases; ACs, adenylate cyclases.
Figure 4 – figure supplement 2 Phylogenetic tree of guanylate cyclase by maximum likelihood (ML). Guanylate cyclase-coupled receptor and soluble guanylate cyclases (sGC) as outgroups. Guanylate cyclases with NIT domains are found in most animal phyla except Porifera, Ctenophora, Urochordata and Chordata. Dot plot of Platynereis NIT-GC genes (columns) expressed in cPRC, INNOS and INRGW (rows) using single cell RNA-Seq. The size of the dots is expressed in proportion to the percentage of cells expressing that gene relative to all cells. The colours represent the normal logarithm of the number of transcripts in the cells expressing the gene.
Figure 4 – figure supplement 3 (A) Co-expression analysis image of the NIT-GC1 (magenta) and MLD-pedal2 amide proneuropeptide gene (MLD: green). Anterior view of the larva at two-day-old. (B, C) Expression analysis of the NIT-GC2 gene (magenta) using in situ HCR. Anterior (B) and posterior (C) views of the larva at three-day-old. (D) Co-localisation analysis using NIT-GC1 (magenta) and NOS (green) antibodies. Anterior view of the larva at two-day-old. (E, F) Localisation analysis using NIT-GC1 and NIT-GC2 antibodies for NIT-GC1 (E) and NIT-GC2 (F) morphant. Green shows co-staining with acetylated α-tubulin antibody (acTub). Anterior view of the larva at two-day-old.
Figure 5 – figure supplement 1 (A) Schematic diagram of the immunostaining procedure after calcium imaging.(B-D) Co-expression analysis image of the NOS (magenta) and RY amide proneuropeptide gene (RYa: green). Anterior view of the larva at two-day-old.